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, We present a variational study of both unpaired and spin-singlet paired states induced in a 

OA ' two-dimensional electron gas at low density by a perpendicular magnetic field. It is based on an 

f-H ^ improved circular-cell approximation which leads to a number of closed analytical results. The 

. ground-state energy of the Wigner crystal containing a single electron per cell in the lowest Landau 

level is obtained as a function of the filling factor v: the results are in good agreement with those 
\^ , of earlier approaches and predict Vc ~ 0.25 for the upper filling factor at which the solid-liquid 

^SJ ■ transition occurs. A novel localized state of spin-singlet electron pairs is examined and found to be 

a competitor of the unpaired state for filling factor v > 1. The corresponding phase boundary is 
quantitatively displayed in the magnetic field-electron density plane. 
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I. INTRODUCTION 



''O ■ The nature of the ground-state of a two-dimensional (2D) many-electron system in a perpendicular magnetic field 
Ch ' has been for a number of years a subject of intense -experimental and theoretical investigation. Since the discovery 
^ , of the Integer and Fractional Quantum Hall EfFectad, this has been one of the richest sources of fundamental new 
physics in Condensed Matter. The transition from a liquid to a Wigner crystal (WC) state, originating from strong 
electron-electron correlations at low density in strong magnetic fields, has received considerable attention (see e.g. 
' ■ Ref. 2 and references given therein). More recently, novel and intriguing behaviors of 2D electron transport in high 
^ \ Landau levels have been reportedcl. In such regime of weak magnetic fields Koulakov et alu predicted that charge 
density waves (CDW) would break the translational symmetry in one direction and form stripes. They also predicted 

■ stabilization of a novel 2D-CDW state with more that one electron per unit cell. Such a "bubble" phase arising from 
> spin-triplet electron pair formation and 2D lattice ordering has been confirmed in systematical numerical studies of 

Haldane et alB and of Shibata and YoshiokaQ. 
OA , In the present work we introduce a new possible low-temperature phase for a 2D electron system in a perpendicular 

■ magnetic field, namely a localized paired-electron state consisting of a spin-singlet electron pair in each cell, and 
4-^ ' investigate its stability at low densities against the more conventional WGnState containing a single electron per cell. 

a . We have in mind situations where the Zeeman splitting can be neglectedQ. In free space, where the Lande (7-factor 
of an electron is 5 = 2, the Zeeman splitting is exactly equal to the cyclotron splitting. It turns out, however, that 
, this is incorrect, for example, in a GaAs heterostructure for the following reasons: (i) the small effective mass in the 
conduction band increases the cyclotron energy by a factor of m/m* ~ 14; (ii) the effective coupling of the electron 
Q spin to the external magnetic field is reduced by spin-orbit scattering by a factor of —5, making the effective Lande 
O ' factor g ~ —0.4. Thus the Zeeman energy is about 70 times smaller than the cyclotron energyQ. 
^ ] This study has two main motivations. Firstly, from the solution of the simple problem of two-electrons moving on 
a plane under parabolic confincmcntEi it emerges that the effective radial potential for the relative motion develops a 
pronounced minimum at sufficiently weak field. The minimum occurs even in the case of vanishing angular momentum 
(to = 0) and is the result of the competition between the Coulomb repulsion and the localizing effects of the magnetic 
field and of the confinement. This feature may hide, in the many-body context and at sufficiently low electron 
densities, some new phase transition to a spin-singlet (for m = 0) paired electron ^state in weak magnetic fields. 
Secondly, and foUowirig, early work on a to = paired state in the 3D electron gasH, such a state was found in a 
variational calculationlld to be a possible competitor to the conventional 3D-WC in zero field. A mixed spin state 
with a preference for spin pairing has been proposed to occur in some regions of the liquid phase from experiments 
on 2D electron systems in the ultra-quantum limittil. 

In our evaluation of the energetics of the paired state we introduce an improved circular-cell approximation for 
an interacting 2D electron gas on a uniform positive background. We first study by this approach the conventional 
unpaired-electron WC state in a magnetic fieldEJ. Each electron is treated as a distributed charge cloud, which can 
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be taken to a good approximation to be localized within a cell if the density is not too high suad the field is not 
too small - the radius of the circular cell being available as an additional variational parameteitJ. Our results are 
in close agreement with those obtained in earlier approaches. We then adopt a similar method to deal with two 
interacting electrons within each circular cell. This yields new and more complex expressions leading to a window of 
thermodynamic stability for the paired-electron state at higher values of the electron density. 

The layout of the paper is briefly as follows. In Section || wc present our treatment of the energy of the unpaired- 
electron 2D-WC in a magnetic field and compare our results with earlier ones and with the energy of the Laughlin 
liquid state. Section III summarizes the main results for the problem of two electrons moving in a plane under an 
applied magnetic field. Section [V reports our treatment of the singlet-paired state and Section ^ gives a variational 
calculation of the energy of this state and a determination of the phase boundary between the unpaired and the singlet- 
paired states. Finally, Section VI summarizes our main conclusions. Some technical points of detail are worked out 
in an Appendix. 



II. DISTRIBUTED-CHARGE APPROACH TO THE WIGNER CRYSTAL IN A MAGNETIC FIELD 

A 2D system of carriers in a pure semiconductor sample subject to a magnetic field B is expected to UB.dergo a 
transition to an ordered triangular lattice structure at sufficiently low temperature T and filling factor itJ. Here, 
1^ = 2-1: UgP where Ug is the areal carrier density and I — {hc/eBy^^ is the magnetic length. The ground-state energy 
of such a 2D-WC has been evaluated by a number of authors with various methods. In the following we shall indicate 
by 6e the values of the ground-state energy after subtraction of the energy of the lowest Landau Level (LO. 

The leading term of a low-density expansion is given by a classical Madelung-potential calculationtj, leading 
to Jeciass = -2.212206 Ryd* = -0.782133 1^^/^ (e^/e^) where e is the dielectric constant of the host medium 
and rsa*g =f—,{Trns)~^^^ with Ryd* and a*^ the effective Rydberg and Bohr radius. Following early Haetree-Fock 
calculationgl3 and a classical-plasma simulation for the determination of the energy-pf the Laughlin liquidtZI, a varia- 
tional calculation based on a correlated magnetophonon wave-function for the WClla gave the result 

5ecwc = -0.782133 0.2410 1/3/2 ^0.16 1/^/2 (e2/e£) (1) 

and predicted that the crystal energy would be lower than that of the liquid only at values of ly lower than ~ 0.14. 
Subsequent work by VignaleuJ used Current-Density Functional Theory in the projected local-density approximation 
and located the upper filling factor for the liquid-solid transition at ~ 0.25, in fairly good agreement with the 
experimental evidence. 

In our approach we model the 2D interacting electron system on a uniform positive background as a collection of 
non-interacting disks, each having a Wigner-Seitz radius rws to be eventually treated as a variational parameter and 
containing a single electron in a trial wave function with a variational width parameter a. Specifically, we adopt a 
Gaussian trial wave function normalized to unity, 

Mr) = (7ra2)-i/2 exp {-ry2a^) (2) 

The single-particle density thus is p{r) ~ <f>sp{r)- Taking the areal density of the background as Pb{r) — UsOir^js — r), 
the total electronic charge inside the disk radius is determined by 

J d^r 4p{r) 0{r„s - r) = 1 - exp {-rlJa') (3) 

and we must require a <C r^s in order to avoid charge leakage errors. As noted by Nagylll, the handling of r^^s as 
a variational parameter ensures that the error involved in setting the electrical potential outside the disk to zero is 
minimized. 



A. Electron-background and background-background interaction energy 



The potential energy Vi„{r) created bsi-the positive background inside the circular cell is given in terms of the Gauss 
hypergeometric function F(a, b; c; d) bytJ 

K„W^^F(l,-i;l;.V.i,) (4) 
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This result is equivalent to the more familiar expression 



(5) 



where E(x) is the complete elliptic integral of the first kind (see Ref. p. 591). 
The self-interaction energy of the background is thus given by 



2 ,/ 37rr„ V r^a 



Ryd\ 



(6) 



having used the result E(x^) xdx = 2/3. Similarly, the interaction energy of the electron with the background is 



ee&(rs,r„s,o-) = 27r / r dr (j)gp{r)Vi,,{r) 
Jo 

^T^sO'B Jo 



For <T ^ r^ivs Eq. (|^) can be replaced by the approximate analytic expression 

" 4 r„s ^ 1 (cr/a|j)2 



dx Ryd* . 



Ryd* 



(7) 



(8) 



This so-called harmonic approximation (HA) for the electron-background energy is obtained from Eq. ^ by extending 
the range of integration up to oo and by using the expansion of the elliptic integral E(a::) ~ (7r/2)(l — x/4) + o{x'^) 
(see|^, p. 591). 

An alternative way to calculate the electron-background energy is by using the expression 



Ceb (^'s , T-ws , 0-) = 27r / pb{r) Vsp (r) r dr , 



(9) 

where Vsp{r) is the electrical potential created by a single electron in the state (||). This can be obtained in closed 
form, 



Vsp{r) = - 



cxp(-rV2a2)Io(rV2(72) Ryd* 



(10) 



where I„(a;) is the modified Bessel function of the ri-th order (see Appendix). The evaluation of Eq. (^ can then be 
carried out in closed form, with the result 
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exp i-rlj2a^) [Mrlj2a') + hirlJ2a')] Ryd* 



(11) 



(see Appendix). The expression (g) is recovered from Eq. ( |ll| ) for rws ^ u by using the asymptotic expansion of 
the Bessel functions (see |2^, p. 377). We shall see below that the "anharmonic corrections" implied by Eq. (11) are 
crucial in a comparison with earlier calculations. 



B. Total energy and variational procedure 



The total energy per electron in a perpendicular magnetic field B is 

et = ek{(T,B) +eeb(rs,rws,o-) -|- ef,b(rs, r„s) 

where tk (f, B) is the kinetic energy. 



(12) 



ek{a,B)^ [ d^r(j)sp{r)-^ (p+ -A)%sp(r) 
J 2m* V c / 



= 27r / r0sp(r) 
Jo 



__l_9.(r5.) + im*a.„V 



(l)sp{r)dr 



(13) 
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with loq = eB/2m*c. We have taken the vector potential in the symmetric gauge, A = B x r/2, and used the fact 
that the wave function is an isotropic state of zero angular momentum. Eq. (^3|) yields 



1 X ^ ^ 

+ As I — 



Ryd* (14) 



\b = m*uj^ay/{2 Ryd*) ~ 3.6 • 10"^^ (em/m*) (B/Tesla)^. For a GaAs heterostructure we have Xb ^ 0.4 
if B = 10 Tesla. 

Equation jl^), after insertion of Eqs. (^), (|l|) and ( |T4| ) is minimized numerically with respect to the variational 
parameters S — cr/a*^ and a = rv^s/Ts<i*B- The equilibrium values of these parameters are shown in Figure |l| as 
functions of Vg for two values of the magnetic field. The ratio "S/tsO. — (cr/r\vs)oq should be appreciably smaller than 
unity for internal self-consistency and the extent to which this consistency criterion is satisfied is shown in Figure |^. 
It is evident from Figures |l| and || that for these values of the field the harmonic approximation is applicable for the 
estimation of the variational parameters whenever the consistency criterion S/r^a <C 1 is satisfied. 

Figure ^ reports our results for the ground-state energy Set as a function of the filling factor in the lowest LL, after 
subtraction of the kinetic energy 2X^J^ Ryd* from the total energy et and rescaling to e^/ei energy units. The left 
panel in Figure ^ shows that Set in these units is still weakly dependent on the field intensity at values of i> larger 
than about 0.4. The right panel in Figure ^compares our results for Set in the range < v < 0.5 with those obtained 
in the correlated Wigner etystal approach of Lam and Girvinllfl (see Eq. (0) ). It is also seen that the classical limit 
of Bonsall and MaradudinEJ is approximately recovered for v ~* 0. 

It is also seen from Figure ^ (left panel) that the ground-state energy of the unpaiced WC as obtained in the present 
approach crosses the energy of the Laughlin liquid as reported by Levesque et alS3 at ly ~ This value for the 

upper critical filling factor of the liquid-solid transition agrees with that reported by Vignalellj and, as discussed by 
this author, there iSnStsong experimental evidence supporting the fact that the crystal exists even at filling factors 
as large as 0.22-0.23a^Eil. However, again as discussed by Vignale, it appears that the liquid-solid transition is not a 
simple crossing of two phases occuring at a single value of v. Rather, it shows a complex reentrant behavior, with the 
liquid phase being stable at or near the Fractional Quantum Hall Effect fractions and the solid phase being stable in 
between. 

We would also like to remark that the simple choice a = I, which ensures that the Wigner-Seitz cell is electrically 
neutral, implies only small changes relative to the prescription proposed by Nagyllj. We in fact find that when we 
take a = 1 the critical filling factor at which the liquid-solid transition occurs shifts from v « 0.25 to ^ 0.23. _ 

Finally, we briefei comment on the issue of Landau Level Mixing (LLM) , as studied in detail by Zhu and LouieEHl 
and by Price et a/Bj. LLM is important when the ratio between the magnetic energy (of order hujo) and the Coulomb 
energy (of order e^/(ersa|j)) becomes smaller then unity, i.e. for ?p-i< A^^. This effect is included in our approach 
through the use of a trial wave function having a variational widthE3 and becomes unimportant as rg increases and 
CTeq approaches the value ^/2£ for which Eq. (g) becomes the exact lowest-energy eigenfunction of the single-particle 
Hamiltonian. This asymptotic behavior is clearly seen from Figure |^ (left panel). 



C. Analytic results in the harmonic approximation 



Having assessed numerically the range of validity of the harmonic approximation (HA) for the estimation of the 
model variational parameters, we proceed to report a number of analytic results which follow from it. The ground-state 
energy is given by 

I]2 16 a^' 

— -I- 



4a 



37rr., 



Ryd* 



(15) 



Minimization of Eq. (^5|) with respect to E and a yields the results shown in Figure ^ 

Let us consider first the choice a — 1 (i.e. = r^a^), where the equilibrium value of S is 



3/4 



(16) 



(l + As r3)i/4 

Thus, Y,a=i{fs, B) decreases with increasing field, due to increased localization of the electron inside the circular cell. 



and saturates to the value A^^^^ corresponding to a^^ = (?i/to*wo)^^'^. The total energy per electron becomes 



^T{rs,B) = 



(1 + XbtIY'^ 



3/2 
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1 



Ryd* 



(17) 
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and increases with the magnetic field, saturating to the value 



B 



2A 



1/2 



Ryd* 



(18) 



The first term in the brackets in Eq. 
vanishing field Eq. yields 



Pq) is the cyclotron zero-point energy hwQ. On the other hand, in the limit of 



tr{r,,B = 0) = 





2 


(stt ^) rs 





Ryd* 



(19) 



2.30 



3/2 



Ryd* 



which is the result obtained by Seidl et aZ@ in their "PC model" (the correlation energy is obtained from Eq. (|l9| ) 
by subtracting from it the exchange energy term given by -8V2/(37rr5) ). The dependence of ef*(rs, S = 0) in Eq. 
( [l9| ) on the density parameter Ts is correct, bttt the numerical coefficients are somewhat different from those which 
are precisely known for the WC in zero fieldll3.|_It may be remarked that the expression for e^^(rs,i? = 0) in 3D 
provides a lower bound to the energy of the WCc3. 

Minimization of the energy in Eq. (|l^) with respect to a is intended to approximately correct for the fact that the 
electrical potential in 2D does not vanish outside the circular cell and is expected to yield a variational lower bound 
for the ground-state energyO. The equilibrium value of the reduced Gaussian width is given by an expression similar 
to Eq. m|). 



Sa(r„B) 



1/4 



(20) 



1 -I- aA_B rf ^ 

while a converges with increasing to the value \/7r/2 (see Figure |l|). That is, Nagy's resultS rws = {V^ I '^)'^ sO'*g 
remains valid in the presence of a magnetic field at low electron density. The expression for the ground-state energy 
becomes 



2%/2 



(l + ABV^r3/2)i/2 4V^i 



rl/4r3/2 



Ryd* 



(21) 



However, this simple analytic expression for the energy of the unpaired WC in a magnetic field is not sufficiently 
accurate on the energy scale needed for the comparisons made in Figure |^. 



III. MOTION OF AN ELECTRON PAIR IN A MAGNETIC FIELD 



In this section we motivate the introduction of a paired state for the WC by discussing the spin-singlet ground-state 
and eflective interaction potential for two electrons moving in a plane under a perpendicular magnetic field. As is 
known frxup the work of Taura, this is an example of a quasi-exactly soluble problem with a hidden s/2-algebraic 
structures: an exact solution exists for special values of the field. Our aim will be to use the analytic form of the 
wave function determined in such a case in order to make a reasonable Ansatz for the variational treatment which 
will be developed for the paired phase in the next Section. 

The Hamiltonian describing the two electrons is written in terms of the relative coordinate r = r2 — ri, of the 
center-of-mass coordinate R = (ri -I- r2)/2 and of the corresponding momenta p and P as 



1 

4m 



-Acm(R) 



1 r 

m 



P + -Ai.ci(r) 

c 



e 
r 



n 



spin 



(22) 



where Acm(R-) = B x R, Ai.ci(r) = B x r/4 and Hspin = —{geh/mc){si + S2) witksi being the spins. The addition of 
a harmonic confining potential merely shifts the value of the cyclotron frequencyQ. The center-of-mass motion in the 
ground-state is described by a Gaussian wave function having width given by the magnetic length £ = {he/ eB)^/'^ . The 
wave function 0roi(r) for the relative motion in the spin-singlet ground-state is even under space inversion. Following 
the treatment developed by TautEl, an analytic expression is obtained for the state of zero relative angular momentum 
(Mrei = 0) if the magnetic field satisfies the condition 7 = 1, where 7 = i-cci/ clb with 4ci — V2^ and aB — fi? /me^. 



This state lies at relative energy 2hjJi^ and apart from a normalization factor is 

0rci(r)oc (l + x)e--'/4 



(23) 
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with X = r/^rei- This is, in fact, the solution which corresponds to the lowest value of 7 and to the lowest energy. 

Two main points of this two-electron problem need emphasizing. Firstly, the effective potential Vcs{i') entering the 
Schrodinger equation for the reduced wave function /eff(r) = <^rei(r) at M^^i = is 



and develops a minimum for 7 > (16/27)1/'' (see Fi gure M). This fact will motivate the introduction of the paired 
phase in the next Section. Secondly, in a solid semiconducting medium the parameter 7 becomes 

7 = — = 7, 25 

a|j em 

7 being the vacuum value. Thus, using — 363.5 (i?/Tesla)~i/^ A and the values e ~ 12.4 and m* /m = 0.067 
for a GaAs heterostructure, we find that the minimum is present in Vcs{r) for all values of magnetic field lower than 
Bcrit - 20 Tesla. 



IV. THE PAIRED PHASE IN A MAGNETIC FIELD 



We return to the circular cell model to develop a variational approach to the total energy in the case of a localized 
paired phase in a spin-singlet configuration. Two electrons are placed inside each disk of radius r^^s , the pair wave 
function being taken in the form 

V'k,<t(i"1, I"2) = Aii.a- (1 + K — —] exp 



4^2 



(26) 



where k and a are variational parameters and AK,a- is a normalization constant. The single-pair case treated in Section 
III is recovered by setting a = I and k = 1/^/2. Notice that the relative motion of the two electrons is described 
in Eq. (|2^) by a Gaussian factor times a linear superposition of the Hermite polynomials of zeroth and first order: 
thus the Ansatz ( p6| ) allows level mixing due to the electron-electron interactions and is the natural extension of 
the variational-width method used for the unpaired phase in Section More refined wave function would include 
higher-order Hermite polynomials in the mixing, weighted by additional variational parameter. 

The main properties of the wave function ( |26| ) are as follows. Firstly, normalization to two electrons per cell yields 



2^1/2 



(27) 



Secondly, the most probable value of the relative distance between the two electrons, as obtained from the appropriate 
maximum of the square of the pair wave function, is 



r2 - ri 



2k 



(28) 



This increases with k and saturates to the value 2a. Thirdly, the single-particle density /9„,cr(r) = J (Pr' \ipf^^„{v,r')\^ 
can be obtained in closed form (see Appendix), with the result 



,-3r^/4(T^ 



2 (^1 + 2^2 + e'''/4-' + « V2^ 



/l->2 ^^2 ;T-.2 ()^2 



^4(7^ 



This expression tends to the correct value Pk=o,(t(^) = (""O"^) ^ exp {—r'^/2a'^) for k — > 0. 

We proceed to evaluate the total energy per electron associated with the wave function (|6|). We have 



(29) 



eb 



^bb 



Impair; I 



(30) 



where for the background-background term we can use the result in Eq. (|g). The electron-background term is 
evaluated, as in the calculation performed in Section A, from the electrical potential VK,a{r) created by the electron 
distribution p^.a-ir) in Eq. (p^) according to 
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eeb^'^TTUs / rdrVf,,cr{r) 
Jo 

4a f°^dx 
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PK,aix) ,h{r„sx/a) Ryd* 



(31) 



(see Appendix), where Jn{x) are Bessel functions. Here 



1 ^ 



1 + 2KV7r + 4k2 



1 + - 



Ke" /8 - 4)Io(a;V8) - x%{x^/8)] 



In the limit of strong locahzation (ct ^ J'ws); we obtain the harmonic- approximation resuh 

8a 2i;2 2 + 5K0F+ 12^2" 



1 + 2Ky/n + 4k2 



(32) 



(33) 



The same hmiting result can be obtained from the electrical potential Kn(^) created by the background disk, which 
is still given by Eq. (||). 

Finally the contribution Cpair in Eq. ( |30| ) is given by 



epair= 2 (V'k,<t|H|V'k,<t) 

,^(^(0rcl|0rcl)(0CM|HcM|0CM) + (0CM|<^CM)('/'rol|Hrel|(?!>rol)) , 



(34) 



the center-of-mass and relative motion states and Hamiltonians being immediately obtained from Eqs. ( p2[ ) and (|2£ 
The result is 



1 2 + 3 K y/TT + 8k'^ 2 2 + 5k y/n + 12k'^ ^/tt I + AK/y/n + 2k' 

' 2Ai3lj ' 



^2I]2 1 + 2kV^ + 4k2 1 + 2kV^ + 4k2 S 1 + 20?/! + 4^2 

The last term in this equation arises from the electron-electron interaction. 



Ryd*. 



(35) 



V. TOTAL ENERGY OF THE PAIRED PHASE 

The total variational energy per electron in the paired circular-cell approximation is obtained from Eqs. (pO|), ( |3] 
( ^ ) and It depends on the three variational parameters S = (j/cig, a = rvi/s/irgO^g) and k. The equilibrium 
values of S and a as functions of show the same trends as displayed in Figure |l| for the unpaired phase. However, 
the asymptotic value of E for an electron pair is reduced by a factor of \/2 and that of a is a = •\/7r/2, increased 
by a factor \/2 over the unpaired state as one expects from the double occupancy of the cell. Internal consistency 
of the theory is ensured by the ratio (cr/rws)oq being smaller than about 0.5, both in the full calculation and in the 
harmonic approximation, for values of > 1 and of the field B > 5 Tesla. The equilibrium value of k as a function 
of Ts is displayed in Figure ^ for two values of the field. The two electrons in each cell are pushed closer together as 
the magnetic field increases, leading from Eq. ( p8| ) to a decrease of k with increasing field as is shown in Figure |[ 

Figure ^ reports the energy Set of the paired phase (that is, after subtraction of the kinetic energy in the lowest 
LL and reduction to e^/(e^) energy units) as a function of the filling factor extending up to v — 1.5. The residual 
dependence of Set on the magnetic field extends into the high-r,, (low-j/) regime, as illustrated in the inset in Figure |^ 
on a magnified scale, and a finite value is attained by Set as ^ 0. These features are a consequence of the Coulomb 
interaction between the two electrons inside each cell. The result in Figure ^ have been fitted to the functional form 
Set — /{el), where 

f{^) = 1/2 ^ f • (36) 

The values of the coefficients in this fitting formula are reported in Table || for various values of the magnetic field. 



7 



Figure [7| compares the energies per electron of the unpaired and paired phases as functions of . The two curves 
cross at two values of Ts for each value of the field, but the crossing at lower is to be discarded as it corresponds 
to situations where charge leakage out of the circular cell is unacceptably high ((cr/rws)cq ~ 1 )• The conclusion from 
the physically significant crossing at larger rj, therefore, is that at each value of the magnetic field in the range from 
5 to 20 Tesla the paired phase becomes stable relative to the unpaired one as decreases. 

From the location of the physically acceptable crossing of the energies of the two phases in Figure we obtain 
their phase boundary in the [vs^Xb) plane, which is reported in Figure |8[ This shows a region of thermodynamic 
stability for the paired phase (at least relative to the unpaired one) in the left portion of the plane, before charge 
leakage outside the cell boundary is expected to occur in an important way so that the model becomes unreliable (to 
the left of the dash-dotted line in the Figure). Withinjiumerical accuracy we find that the boundary between the 
two phases is actually set by the v = 1 line in the planeE3. Therefore, the present simple model predicts that electron 
clusterization may occur only in Landau levels above the lowest one. Of course, the unpaired state in the v < 1 region 
(on the right of the v — 1 line) could become stable against the Laughlin liquid only for v < 0.25, as already discussed 
in connection with Figure |[ 

We should emphasize that the ground-state energy difference that we estimate between the paired and the unpaired 
state is quite appreciable as long as the magnetic field intensity lies in the range that we have illustrated in our 
calculations. The difference decreases with the field intensity, so that more refined calculations would be needed at 
low fields. 



We have presented a variational study which approximately accounts for the ground-state energetics of a 2D many- 
electron system with Coulomb interactions at low density and in the presence of a perpendicular magnetic field. The 
method consists of an improved 2D version of the well-known spherical cell approximation and is applied to two 
distinct situations, a single-electron state and a spin-singlet paired state. For the former it gives new analytical and 
numerical results which are in close agreement with the state-of-the-art calculations on the 2D Wigner crystal. 

The evaluation of the spin-paired state has been motivated by previous studies and suggests that this state may be 
stable in a range of system parameters corresponding to Landau levels above the lowest one. It remains to be seen 
whether such a paired state would be confirmed in more sophisticated treatments allowing in particular for inter-cell 
corrections. It may be also emphasized at this point that we have not examined the stability of the spin-singlet paired 
state against thepemergence of spin-polarized states. Evidence from measurements of Knight shift of the ''^Ga nuclei 
in n-doped GaAsE3 indicates that this quantity, which is proportional to the spin polarization, drops precipitously on 
either side of = 1, which is evidence that the charged excitations of the 1^=1 ground state are finite-size Skyrmions. 
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APPENDIX. DETAILS OF ANALYTIC CALCULATIONS 



We report in this Appendix some details on the derivation of some analytic results given in the main text. 



A. Equations and (|ll|) 



The single-electron potential Vsp (r) is 




(37) 



(in Ryd*), where 




(38) 
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Equation ( p^ in the main text follows by using the result 



exp (-a;74) Jo(a x) ^ ^/^ cxp (-072) Io(a72) , 



Finally, the integral in Eq. @ is carried out with the help of the relation 



X exp (-.tV2) dx^'- er'"/^ [UbV2) + Ii(6V2)] 



leading to Eq. (nil) in the main text. 



(39) 



(40) 



B. Equation (29) 



The single-particle density pK^ui^) is calculated from the expression 

Pk,o-(?') = 2f7^^^ exp (-r^/cr^) / d0 x dx {1 + kx)"^ e^p {-x r cos9/a - /2) . 

Jo Jo 



and 



Eq. (pol ) is obtained by using the results 



/ exp (— xz/cos^) d0 = 27rIo(x2/) 
Jo 



+00 





exp 



X (1 + Kx)^ exp (— a;^/2) Io(a;?;) dx = 
(2/V4) { (1 + + kY) exp (yV4) + [(2 + 2/') v^Io(yV4) + 2/' Ii(yV4)] } 



(41) 



(42) 



(43) 



C. Equation (|^) 

The electrical potential VK.air) created by the electron distribution pK.air) in Eq. (E9h is written as 



+00 



(44) 



(in Ryd*), where in essence pK,aix) is the Fourier transform of the electron density distribution (see Eq. ( ^ ) ). After 
inserting this formula into the first line of Eq. ( |3l| ) and interchanging the order of the two integrations, the integration 
over r can be carried out with the help of the result 



' b 
Jq {ax) X dx ^ — J i{a b) . 



(45) 



This yields the second line of Eq. IpW) in the main text 
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TABLE I. Coefficients of the interpolation formula in Eq. ( 



B (Tesla) 


a 


b 


c 


d 


e 


/ 


5 


-0.6621 


0.2769 


-0.0647 


0.1751 


-0.0072 


0.5578 


10 


-0.6435 


0.2715 


-0.0588 


0.1828 


-0.0231 


0.5406 


15 


-0.6325 


0.2655 


-0.0548 


0.1869 


-0.0288 


0.5302 


20 


-0.6255 


0.2612 


-0.0522 


0.1894 


-0.0315 


0.5237 
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FIG. 1. Equilibrium values of the reduced width E = (a |ag)^^^ (left panel) and of the reduced disk radius a = (rws/''s<iB)cq 
(right panel) as functions of Vs and for two values of the massjetic field. The asymptotic values at large are acq = (ft/m'ojo)^''^ 
and a = the latter being the value obtained by Nagytj in zero field. The dot-dashed lines are obtained from the harmonic 

approximation, which is becoming approximately valid for > 1. 
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FIG. 3. Left panel: Ground-state energy Set of the unpaired WC, in units of e^/e€, as a function of the filHng factor v in 
the range < i/ < 1 for values of B = 5, 10, 15 and 20 Tesla (from bottom to top). Right panel: Ground-state energy of the 
unpaired WC obtained in the present approach (full line) compared with the CWC result of Lam and Girvin (dotted line), 
with the asymptotic classical result (dashed line) and with the energy of the Laughlin liquid (dot-dashed line, from Levesque 
et al.). 
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FIG. 4. Effective potential 14ff(r) for the relative motion at zero relative angular momentum. The dotted line corre- 
sponds to the critical value 7crit = (16/27)^''*, above which a minimum is present. The dotted line has an inflection point at 

Xciit = 2/(37crit). 



15 




FIG. 5. The parameter k as a function of rs for two values of the field B in the paired phase (full lines). The dash-dotted 
lines show the results of the harmonic approximation. 



16 




17 




FIG. 7. Comparison of the total energies of the paired state (full lines) and of the unpaired state (dash-dotted lines) for two 
values of the magnetic field. 
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FIG. 8. Boundary between the singlet-spin paired state and the unpaired state in the (rs,As) plane. The boundary is set 
by the z/ = 1 line. The dash-dotted line shows where the consistency ratio E/rsO; for the paired state becomes unity. 
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